Rapidly rotating spherical Couette flow in a 
dipolar magnetic field: an experimental study 
of the mean axisymmetric flow 



Nataf, H.-C. * Alboussiere, T., Brito, D., Cardin, P., 
Gagniere, N., Jault, D., and Schmitt, D. 

Geodynamo team, LGIT-UMR5559-CNRS-UJF, Grenoble, France 



Abstract 

In order to explore the magnetostrophic regime expected for planetary cores, in 
which the Lorentz forces balance the Coriolis forces, experiments have been con- 
ducted in a rotating sphere filled with liquid sodium, with an imposed dipolar mag- 
netic field (the DTS setup). The field is produced by a permanent magnet enclosed 
in an inner sphere, which can rotate at a separate rate, producing a spherical Cou- 
ette flow. The flow properties are investigated by measuring electric potentials on 
the outer sphere, the induced magnetic field in the laboratory frame just above the 
rotating outer sphere, and velocity profiles inside the liquid sodium using ultrasonic 
Doppler velocimetry. The present article focuses on the time-averaged axisymmetric 
part of the flow. The electric potential differences measured at several latitudes can 
be linked to azimuthal velocities, and are indeed found to be proportional to the 
azimuthal velocities measured by Doppler velocimetry. The Doppler profiles show 
that the angular velocity of the fluid is relatively uniform in most of the fluid shell, 
but rises near the inner sphere, revealing the presence of a "magnetic wind", and 
gently drops towards the outer sphere. The transition from a magnetostrophic flow 
near the inner sphere to a geostrophic flow near the outer sphere is controlled by 
the local Elsasser number. For Rossby numbers up to order 1, the observed velocity 
profiles all show a similar shape. Numerical simulations in the linear regime are 
computed, and synthetic velocity profiles are compared with the measured ones. A 
good agreement is found for the angular velocity profiles. In the geostrophic region, 
a torque-balance model provides very good predictions. Radial velocities change 
sign with the Rossby number, as expected for an Ekman-pumping dominated flow. 
For a given Rossby number the amplitude of the measured angular velocity is found 
to vary by as much as a factor of 3. Comparison with numerical simulations sug- 
gests that this is due to variations in the electric coupling between liquid sodium 
and the inner copper sphere, implying an effect equivalent to a reduction of the 
inner sphere electric conductivity by as much as a factor 100. We show that the 
measured electric potential difference can be used as a proxy of the actual fluid ve- 
locity. Using this proxy in place of the imposed differential velocity, we find that the 
induced magnetic field varies in a consistent fashion, and displays a peculiar peak 
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in the counter-rotating regime. This happens when the fluid rotation rate is almost 
equal and opposite to the outer sphere rotation rate. The fluid is then almost at rest 
in the laboratory frame, and the Proudman-Taylor constraint vanishes, enabling a 
strong meridional flow. We suggest that dynamo action might be favored in such a 
situation. 

Key words: Spherical Couette flow, magnetostrophic, liquid sodium experiment, 
dynamo, Taylor-state, ultrasonic Doppler velocimetry. 
PACS: 47, 47.65, 91, 91.25 



1 Introduction 



The internal dynamics of liquid planetary cores is deeply influenced by the 
rotation of the planet and (in many cases) by the presence of a magnetic 
field, and can be in the so-called magnetostrophic regime where the Coriolis 
and Lorentz forces balance each other. In recent years, attention has been 
paid to the rotating spherical Couette flow in the presence of a magnetic 



field (iKleeorin et all. 119971; IHollerbachl. Il994l ; iDormy et all . Il998l ; IStarchenko 



19981 ; IHollerbachl . l2000l ; iNataf et all 120061 ) . An electrically conducting liquid 



fills the gap between two concentric spheres and is sheared through the dif- 
ferential rotation of the two spheres. A magnetic field can be applied. This 
flow has the distinct advantage of being amenable to both experiments and 
numerical models, whil e retaining important ingredien t s of n atural situations. 
Numerical simu l ation s ( Dormv et al.,ll998"l; Hollerbach . 2000 ) and asymptotics 
f lKleeorin et all Il997l : IStarchenkol . ll998l : iDormv et all l2002h highlight the im- 
portance of the magnetic wind that arises when the magnetic field lines are 
not parallel to the rotation cLXIS, clS for a dipolar field. This magnetic wind 
can entrain the fluid at angular velocit ies larger than that o f the inner sphere, 
a phenomenon called super-rotation (IDormy et all Il998l . 120021 ). The rotat- 
i ng spherical Couet te flow has also been proposed for producing a dynamo 
( jCardin et all 120021 ). Numerical models have shown that dynamo action could 
be produced this way , even in liquids with a low magnetic Prandtl number, 
such as liquid metals (jSchaeffer and Cardinl . 120061 ). This has led Dan Lathrop 
and his team to build a 3m-diameter sphere with a rotating inner sphere in 
the hope of starting a dynamo, once the gap is filled with liquid sodium and 
the differential rotation is large enough. 
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Here, we report on results obtained in th e DTS (Derviche To urn eur Sodium) 



experi ment. This experiment, described in lCardin et al.l (120021 ) and lNataf et al. 
( 120061 ). is a rotating spherical Couette flow with an imposed axisymmetric 



dipolar magnetic field. Because of this strong imposed field, and although 
liquid sodium is used, the dynamo effect cannot be studied in the DTS ex- 
periment. However, it permits to explore very similar dynamics. In particular, 
the magnetic field is large enough to strongly influence the flow (interaction 
parameter N up to 200) and the magnetic field is strongly affected by the flow 
(magnetic Reynolds number Rm up to 44). 



This kind of experiment can bring information on the level and organiza- 
tion of turbulence in situations where both rotation and the magnetic field 
i mpose strong constraints on the flow. Local theories have been pr oposed 
( Braginsky and Meytlie . 1990h . partly backed by numerical simulations ( St. Pierre 
19961 ). and some teams are tryin g to imp l emen t sub-grid models based on these 
theories in full dynamo models (IBuffettl . 120031 ; iMatsui and Buffettl . 120071 ). Re- 
sults obtained in the D TS experiment suggest an alternative view, in which 
waves play a major role (ISchmitt et al.l . 120081 ). A good characterization of the 
mean axisymmetric flow is needed in order to fully exploit these observations. 



An experimental se tup, built by Dan Lathrop and his group (ISisan et al.l . l2004l ; 
Kelley et al.l . 120071 ) is very similar to our DTS experiment. The main dif f erence 
is th at an axial magneti c field is applied by external coils in lSisan et al.l (120041 ) 
and iKelley et al.l (120071 ). while we impose a strong dipolar magnetic field in 
DTS since the inner sphere contains a permanent magnet. One of the conse- 
quences is that the entrainment of the conducting fluid is much stronger, since 
the inner sphere acts as a magnetic stirrer. We also have developed specific 
instrumentation, such as measuring the electric potential at the s urface of the 



outer sphere, which carrie s important information on the flow (INataf et al. 



20061 : ISchmitt et all I2008T ). 



In this article, we focus on the time-averaged axisymmetric magnetohydrody- 
namic flow observed in DTS when the oute r sphere rotates. The first obser- 
vations were described in lNataf et al.l (120061 ). and more thorough results for a 
static outer sphere will be discussed in a forthcoming article. 



Nataf et al.l (120061 ) showed that the measurement of electric potential differ- 
ences at the surface of the outer sphere could be related to the azimuthal 
velocity of the liquid beyond the Hartmann boundary layer. It was found that 
this velocity could be larger than the solid-body rotation of the inner sphere, 
pr oviding exper i menta l evidence for the so-called super-rotation predicted 
by iDormy et al.l (119981 ). However, the latitudinal variation of the electric po- 
tential differences was found to differ marked ly from the o n e calc ulated in a 
linear numerical model, analogous to that of IDormy et al.l (119981 ). The rela- 
tive variation of electric potential differences with latitude turned out to be 
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universal, independent of the differential rotation of the inner sphere. This 
was particularly true when the outer sphere was at rest. The electric design of 
the electrodes and acquisition system has been improved since, yielding better 
quality signals. 



Numerical models provide a useful guide for interpr eting; the expe r iment al ob- 
servations. Linear solutions have been computed by iDormy et al.l (119981 ) for a 
situation similar to ours (spherical Couette with an imposed dipolar magnetic 
field, a conducting inner sphere and a rapidly rotating outer sphere). The re- 
sults show that as the magn etic field is i ncrea sed, the flow evolves from the 
classical Stewartson solution (IStewartsonl . Il966l ) . with a cylindrical shear-layer 
attached to the rotating inner sphere, to a solution where the fluid enclosed 
within a closed magnetic field line is entrained by the inner sphere to angular 
velocities in excess of the imposed differentia l rotation. A shear-layer devel- 
ops around that field line (IDormy et al.l . |2002| ). Using the numerical approach 
described in Nataf et al. (l2006h . we have extended these computations to con- 
ditions closer to the actual DTS configuration, with stronger global (/) and 
differential (A/) rotation rates. The results are presented in figure [TJ The 
angular velocity contours (figure dh) clearly illustrate the existence of two do- 
mains: a region close to the inner sphere where the magnetic effects dominate, 
and a region where the flow is geostrophic farther away. In the inner region, 
the angular velocity conto urs follow the magnetic field lines (Ferraro law of 
isorotation (jFerrard . I1937T) ^ . while they ar e aligned with the rotation axis in 
the geostrophic region. IDormy et al.l (119981 ) demonstrate that the radial extent 
of the magnetic-dominated region solely depends upon the Elsasser number 
(which compares the Lorentz and Coriolis forces). In contrast to the Stew- 
artson solution, the bulk of the fluid is entrained by the inner sphere, and a 
thin Ekman-Hartmann layer forms beneath the outer shell. Ekman pumping 
results and feeds the meridional circulation shown in figure [Tb. The circula- 
tion around the main vortex is counterclockwise in the upper hemisphere [i.e. 
centrifugal at the equator) when the inner and outer spheres rotate in the 
same direction, and clockwise {i.e. centripetal at the equator) when the two 
spheres are in contra-rotation. The field lines of the meridional magnetic field 
produced by this circulation are drawn in figure [lb. 



Direct measurements of the flow velocity were presented in lNataf et al.l (120061 ). 
using ultrasonic Doppler velocimetry, but only radial profiles of the radial ve- 
locity could be obtained. More profiles are now available and can be compared 
with the numerical simulations. The setup now permits the record of profiles 
of the angular velocity as well, revealing the actual variation of angular ve- 
locity inside the sphere. We will see that they confirm the progressive rise of 
the angular velocity from the equator of the outer sphere to the interior of 
the fluid predicted by the numerical simulatio ns (figure [Th^ . The ex istence of 
this peculiar boundary layer was predicted by iKleeorin et al.l (119971 ). We will 
present their approach and propose an extension that provides a better fit to 



4 



our observations. 



There were some hints in iNataf et al.l (120061 ) that different flows could be 



observed for the same imposed velocities. With additional data now available, 
this tendency is fully confirmed. We think that this is due to variations in 
the electric coupling between the liquid sodium and the copper inner sphere. 
Indeed, the entrainment of the liquid is not viscous but depends on electric 
currents that flow between the liquid sodium and the copper shell of the inner 
sphere. The presence of oxides or other impurities, or imperfect wetting on 
that shell can hinder the electric currents. However, we show in this article 
that we can use the difference in electric potentials at one latitude (40° for 
example) as a proxy of the actual differential angular velocity of the fluid, 
and get a much more coherent picture of the various observations, if we relate 
them to this actual value. 



No m easurement of the induced magnetic field was provided by INataf et al. 



(120061 ). Here, we present measurements of the radial and orthoradial compo- 
nents of the induced magnetic field at a latitude of about 50° just outside the 
outer sphere, in the laboratory frame. The mean values correspond to the field 
induced by the mean meridional flow. When the outer sphere is rotating, a 
very peculiar behavior is monitored for a particular counter-rotation rate of 
the inner sphere. 

We recall the main features of the DTS setup and measuring techniques in 
section [2j The typical dimensionless numbers are discussed in section [3j The 
results are presented in section 0], and discussed in section 



2 Experimental setup and measuring techniques 



2.1 Setup 



The central part of the DTS experiment is sketched in figure [2J Forty liters of 
liquid sodium are contained between an inner sphere of radius b = 7.4 cm and 
an outer spherical surface of radius a = 21 cm. The 5 mm-thick outer sphere 
is made of stainless steel, which is 8 times less electrically conductive than 
sodium at 130°C, the typical experimental temperature. The inner sphere is 
made of two hollow copper hemispheres filled with 5 layers of magnetized rare- 
earth cobalt bricks. Copper is about 4.2 times more conductive than sodium. 
Outside the copper sphere, the magnetic field is found to be an axial dipole, 
aligned with the axis of rotation and pointing downwards with: 

B(r,6) = ^^(2e r cosfl + e e sin#), 
47rr d 
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Fig. 1. Non-linear axisymmetric numerical simulation with the geometry and pa- 
rameters of the DTS setup, corresponding to / = 5 Hz and A/ = 0.4 Hz (see section 
[3] below) (a) Contour map of the angular velocity. The interval between contours 
is 0.027 in 2irAf units. Angular velocities are larger than that of the inner sphere 
in the shaded region (super-rotation). Two different regions are clearly visible: the 
region near the inner sphere, where the flow lines are aligned with the magnetic 
field lines of the imposed dipole (Ferraro law), and where the angular velocity ex- 
ceeds that of the inner sphere by up to 6%; farther away, a region where the flow 
is nearly geostrophic, with contours aligned with the rotation axis (vertical). Also 
drawn, the trajectory of the ultrasonic beam shot from the "azimuthal" assembly 
at +10°, projected on the meridional plane. Marks are drawn every 2 cm along the 
ray path, (b) Meridional circulation. The interval between contours is 2.7 x 10~ 5 
in 27ra 2 A/ units, where a is the radius of the outer sphere. The main circulation 
is counterclockwise, i.e. centrifugal at the equator. It changes sign with A/. Radial 
velocities are measured along radial shots at latitudes 40° and 20° (lines), (c) Field 
lines of the induced meridional magnetic field. The interval between contours is 
2 x 10~ 3 in aBo units. Bq is the intensity of the imposed dipolar magnetic field at 
the equator at the outer sphere. The line is thicker. While the imposed dipole is 
pointing downwards, the induced field points up at high latitudes. B r and Bq are 
measured at a latitude of 50°, 1 cm away from the sodium boundary (arrows). 
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in spherical coordinates, with 9 the colatitude and Ai = —700 Am 2 the mag- 
netic dipolar moment, yielding a magnetic field amplitude ranging from 0.345 
T at the poles of the inner sphere down to 0.008 T a t the equator of the outer 
sphere. More details are given in iNataf et al.l (120061 ). 



The inner and outer spheres are set in rotation around the vertical axis. In 
this article, we present results obtained for imposed rotation rates / of the 
outer sphere of 5, 10 and 15 Hz. The rotation rate of the inner sphere in 
the laboratory frame is varied between -25 and +25 Hz, yielding differential 
rotation rates Af of the inner sphere in the rotating frame ranging from -40 
to 20 Hz (depending upon /). Each sphere is entrained by an 11 kW brushless 
electric motor (Yaskawa SGMGH-1ADCA61). 

The torques and actual rotation rates are retrieved as analog signals from the 
two motor controllers. In order to decipher the organisation of the flow in the 
liquid, several techniques are employed. Electric potentials are measured at 
several locations on the outer sphere. Profiles of the flow velocity are obtained 
using ultrasonic Doppler velocimetry. Since the outer sphere is rotating, all 
these signals must pass through electric slip-rings to be analyzed in the labo- 
ratory frame. We use a 36-tracks 10A slip-ring system (Litton AC6275-12). 
In order to decrease electromagnetic perturbations of the signals, this sys- 
tem is completely enclosed in an aluminum casing. In contrast, the magnetic 
field is measured in the laboratory frame, with magnetometer sensors a few 
millimeters away from the rotating surface of the outer sphere. 



2.2 Electric potentials 



We measure the electric potential at several points on the outer sphere. Blind 
holes are drilled in the 5 mm-thick stainless steel shell. The holes are 1 mm in 
diameter, and 4 mm-deep. They are threaded and equipped with 3 mm-long 
Ml. 2 mm brass bolts, on the head of which electrodes are soldered. Stainless 
steel is a poor electric conductor as compared to liquid sodium and thus it 
does not affect the dynamically generated electric potential, but it is a good 
enough electric conductor for the impedance of electrode pairs to be much 
smaller than that of the acquisition system. The measurements are acquired 
at 1 kHz with a PXI-6229 National Instruments board, after going through a 
simple anti-aliasing 215 Hz low-pass RC filter. As we focus here on the mean 
axisymmetric state, these data are time-averaged over time-windows ranging 
from 0.1 to 1 second, in order to remove fluctuations. 

In this article, we present data for the difference in electric potential measured 
between electrodes placed along the same meridian. These electrodes are 10° 
apart and range from —45° to +45° in latitude as shown in figure [21 Assuming 
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Fig. 2. Sketch of the central part of the DTS experiment. The 7.4 cm-radius inner 
sphere is made of copper and contains 5 layers of magnets. The inner radius of 
the stainless steel outer sphere is 21 cm and its thickness is 5 mm. Blind holes are 
drilled and threaded along several meridians and parallels to receive Ml. 2 mm brass 
bolts acting as electrodes. In this article, results are presented for electric potentials 
measured along one meridian (on the left). The large holes at latitudes —20°, 10° 
and ±40° receive interchangeable assemblies, which can be equipped with ultrasonic 
transducers. Magnetometers are placed at a latitude of 50°, just above the outer 
sphere, in the laboratory frame. They measure the radial B r and orthoradial Bg 
components of the magnetic field. 



steady state, the difference in electric potential AV^ between two electrodes 
separated by a latitudinal angle AO is related to the azimuthal velocity U<p of 
the liquid sodium by: 

AV „ „ jg 



U^Br 



T 



aA9 a 

where B r is the radial component of the magnetic field, jg is the ^-component 
of the electric current, and a is the electric conductivity of liquid sodium. 



In iNataf et al.l (120061 ). we used the differences in electric potential to deduce 
the angular velocity of the liquid sodium flow beneath the Hartmann boundary 
layer. Indeed, assuming that the meridional electric currents vanish below that 
layer, equation [1] yields: 



AV 
aA9 B r 



or A/ fluid (s) 



AV 



2ns a AO B r 



(2) 
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where A/fl ui d(s) is the angular frequency of the fluid with respect to the outer 
sphere at the cylindrical radius s for the latitude of the electrode pair. In this 
article, we will use the difference in electric potential measured at 40° latitude 
(AV4 °) as a proxy of the actual rotation frequency of the liquid sodium flow 
with respect to the outer sphere. 



2.3 Ultrasonic Doppler velocimetry 



Since liquid sodium is opaque, it is not possible to use optical methods to inves- 
tigate the flow within the sphere. Instead, we use ultrasounds to probe the flow. 
Ultrasonic transducers (TR0405HS from Signal Processing, Switzerland) are 
placed in removable assemblies on the outer sphere, and shoot a narrow beam 
of 4 MHz ultrasounds. The pulsed Doppler velocimetry technique is perfectly 
suited for retrieving profiles of the component o f the flow velocity al o ng; th e 
sho oting line, as illustr ated in lBrito et al.l (120011 ) . lEckert and Gerbethl (120021 ). 
and iNataf et al.l (120061 ). We use a DOP 2000 instrument from Signal Process- 
ing, Switzerland. 



While only radial profiles of the radial velocity were discussed in INataf et al. 



(120061 ). we show here profiles of both the radial and angular velocities. Radial 
profiles are measured at two different latitudes (—20° and —40°), as drawn 
in figure [lb. The angular velocities are obtained with an assembly devised 
for shooting a beam at an angle of 24° away from radial. Radial velocities 
are found to be at least one order of magnitude smaller than the azimuthal 
velocities. Neglecting the radial velocity component, the axisymmetric angular 
frequency A/ fluid (s, z) of the fluid flow can thus be retrieved by projection of 
the velocity U (d) measured along the shooting line. With our actual geometry, 
one gets: 



A,/fl ui d 0, z) 



2.88W(d) 
2tt a 



(3) 



where d is the distance along the shooting line, s the cylindrical radius, and z 
the height above the equatorial plane, which are related by: 



z(d) = a sin 10° - 0.359 d 
s(d) = a cos 10° 



\ 



d I d 

1 - 1.76- + 0.898 - 
a \ a . 



(4) 
(5) 



As shown in figure [Th^ plotting z versus s for this beam yields the projection 
of the ray path in a meridional plane. Tick marks have been drawn every 2 
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cm along the path. The minimum cylindrical radius reached is s = 7.7 cm, 
close to the radius of the inner sphere (7.4 cm), but it is attained some 4 cm 
beneath the equator. 

The Doppler technique relies upon ultrasonic energy backscattered by particles 
in the liquid. We did not add any particle in the liquid sodium. In liquid metals, 
it is often assumed that the backscattered echoes are due to floating oxides. 
It is interesting to note that we were able to record good profiles, even when 
the outer sphere was rotating at 14 Hz, implying centrifugal acceleration up 
to 180 times the acceleration due to gravity. 

2.4 Induced magnetic field 

The magnetic field is measured outside the sphere, in the laboratory frame. 
Both the radial B r and the orthoradial Bq components are measured with a 
"Giant MagnetoResistance" (GMR) chip. The two chips are mounted on a 
circuit protected by a stainless steel tube. The probe has been installed at 
a latitude of 50°, as close as possible to the rotating sphere: the distance to 
the center of the sphere is about 22 cm. The variation of temperature has a 
significant effect on the measurements. Temperature is thus measured close 
to the GMR magnetometers, and a linear correction scheme is applied. The 
acquisition rate is 2000 samples per second, but all data presented here have 
been moving-averaged over 3 seconds, since we are focusing on the time- 
averaged axisymmetric mean flow. 



3 Dimensionless numbers 

Relevant parameters and dimensionless numbers are given in table [TJ The 
symbols have their usual meaning : rj is the magnetic diffusivity, p the density, 
v the kinematic viscosity of the liquid sodium. As in all moderate-size liquid 
sodium experiments, the Joule dissipation time tj is relatively short. Since the 
imposed magnetic field is dipolar, its amplitude is much larger near the inner 
sphere (B = Bj) than near the outer sphere (B = B ). We therefore evaluate 
those numbers that depend on the strength of the magnetic field in these 
two regions. The Hartmann number Ha is large everywhere. The Lundquist 
number S is large near the inner sphere, implying that Alfven waves can play 
a role in the dynamics, with typical velocities U a ranging from 0.2 to 5 m/s, 
in the same range as flow velocities. 

Values of the Ekman number E are given for two typical values of the rota- 
tion rate / of the outer sphere. The DTS experiment has been designed for 
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the study of the magnetostrophic regime, in which Lorentz forces and Coriolis 
forces are of similar amplitude. This is measured by the Elsasser number A, 
which ranges from 0.01 near the outer sphere (where Coriolis forces domi- 
nate) to nearly 10 near the inner sphere (where Lorentz forces dominate). For 
time-scales short compared to the magnetic diffusion time, a measure of the 
relative influence of the rotational and magnetic effects is given by the dimen- 



sionless number A, t 



Fearn et al.. 


1988; 


' Lehnert . 


195^ 


1) bv 



re ratio of the Alfven velocity to the inertial wave velocity 
Cardin et all |2002| ). which is called the Lehnert number 



Jaultl (120081 ). The Lehnert number is small in planetary 
cores, and rotational effects are expected to dominate at short time-scales. In 
the DTS experiment, the Lehnert number is less than 1. 

We estimate a typical velocity U from the magnitude of the differential rota- 
tion Af of the inner sphere with respect to the outer sphere. The estimate is 
based on the tangential velocity at the equator of the inner sphere. We will 
see that, because the magnetic coupling between the magnetized inner sphere 
and the liquid sodium is strong, velocities in the fluid layer typically reach this 
value. This yields magnetic Reynolds numbers Rm in excess of 40, while the 
Reynolds number Re lies in the range 10 5 — 10 6 . The interaction parameter N 
ranges from 0.01 near the outer sphere to 250 near the inner sphere, meaning 
that inertial effects only play a role near the outer boundary. Note that the 
Rossby number Ro is here simply defined as Ro = AQ/Q = Af / f . 

The DTS experiment thus enables us to explore an original parameter range, 
which covers some aspects of liquid core dynamics. 



4 Results 



4-1 A typical run 



Figure [3] illustrates the data typically acquired during a 1600s-long run. The 
top graph shows the variation of /, the rotation frequency of the outer sphere, 
and Af, that of the inner sphere with respect to the outer sphere. The outer 
sphere is first set to a rotation rate of about 5 Hz, then several plateaux 
of different Af are produced, down to Af = -21.5 Hz. Here, they all corre- 
spond to counter-rotation. Some small fluctuations of the rotation rates are 
observed (±3% max for Af and ±7% max for /). They are due to variations 
in the power demanded to the motors, as shown in the next graph below, 
which displays the total power of the motors, reaching 2.3 kW. Note the large 
power fluctuations (±22%) in the middle of the plateau with Af = -19 Hz. 
The next graph gives the differences in electric potential AV recorded by the 
eight electrode pairs along a meridian. The highest signal reaches 3 mV. As 
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Fig. 3. An example of the data acquisition during an experimental run, as a function 
of time (in seconds). The top graph shows /, the rotation rate of the outer sphere, 
and A/, the rotation rate of the inner sphere with respect to the outer sphere. The 
next graph below is the total power. The differences in electric potential recorded 
by the eight longitudinal electrode pairs (from —40° to 40° latitude) is shown next. 
The last graph below gives the B r and Bg components of the induced magnetic field 
(in mT) recorded on the 50°datitude magnetometer, corrected for temperature. The 
horizontal bar in that graph gives the time window of the angular velocity record 
shown in figured! 
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symbol 


expression 


units 




value 


a 


outer radius 


cm 




21 





lllilcl I dU.1 Uo 


cm 




7 A 


Tj 


/--, 2 / 2 /v-, 
Ct / 71 Tj 


s 


B = Bi 


u.uo 

B = Bq 


Ha 


aB/^fiopvr] 




4400 


210 


Q 

o 






1 


U.OD 


TT 


r / rr, — 7i 
B/^/PoP 


m s~ 


5.1 


0.2 








/ = 5Hz 


/ = 15Hz 


T ' 

hj 


, . /r>„2 




4.7 10" 7 


1.6 10" 7 








5 = 5, S = B B = Bi B = B 


A 

A 


_ rj2 / _r~> 

cri3 / ps I 




9.2 0.02 


3.1 0.01 


A 






0.77 0.04 


0.26 0.01 








Af = 0.5Hz 


Af = 40Hz 


J T 

u 




m s 


0.23 


18.6 


xim 


Ua/rj 




0.55 


44 


Re 


Ua/v 




7 10 4 


6 10 6 








i? = Bi B = B B = Bi B = B 


N 


aaB 2 / P U 




262 0.06 


3.1 0.01 



Table 1 

Typical values of the relevant parameters and dimensionless numbers for different 
imposed rotation frequencies / = 0,/2-k of the outer sphere and differential rotation 
frequencies Af = A£1/2tt of the inner sphere. For the numbers that depend on the 
magnetic field strength, two values are given, the first one with B = Bi = 0.175 T 
at the equator of the inner sphere, the second one with B = B Q = 0.008 T at the 
equator of the outer sphere. 

expected from equation [2], the potential difference is largest at the highest lat- 
itude because B r is largest there. Since Af is negative, the electric potential 
differences are positive in the southern hemisphere and negative in the north- 
ern hemisphere (with our connection conventions). There is an almost perfect 
symmetry between the two hemispheres. Several ultrasonic Doppler velocity 
profiles are recorded during this sequence. The last graph below shows the 
induced magnetic field measured at 50° latitude. The records have been time- 
averaged along a 3 seconds-window to remove the fluctuations linked with the 
rotation of the outer sphere, and the imposed dipole contribution has been re- 
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Fig. 4. An example of an ultrasonic Doppler angular velocity acquisition. The hor- 
izontal axis is time (in seconds) and the vertical axis is distance along the shooting 
line (in cm), here down to d ~ 10 cm. The corresponding time-window is drawn in 
figure [3l The shades give the velocity value in m/s (scale on the right-hand side). 
Strong velocity fluctuations are observed. 



moved. The plateaux are well identified, but large fluctuations are observed 
when Af is large. The maximum amplitude of the induced magnetic field is 
0.23 mT (2.4% of the imposed field at this position). 



A typical ultrasonic Doppler velocity record is shown in figure HI It corresponds 
to the time window marked with a horizontal bar in the last graph of figure 
|3l when time-oscillations are observed on all data. The spatiotemporal plot 
clearly shows these oscillations. The velocities are largest at the bottom of 
the figure, where the ultrasonic beam gets closest to the inner sphere for this 
record. Equation Ogives the value of the cylindrical radius s = 12.1 cm. There, 
the velocity reaches -3.5 m/s. This large value is nevertheless smaller than the 
imposed tangential velocity on the inner sphere (-8.7 m/s). For this file, the 
time-interval between profiles is 15 ms, and there are 128 shots per profile. 
Velocities cannot be measured in the first 1 cm of the records, due to strong 
reverberations in the assembly wall. 



Note that all data display large fluctuations in the time- window between 1100 
and 1250 seconds. We will see that, in this parameter range, the flow exhibits a 
fundamental transition between two different configurations (see section l5~4j) . 
Here, this transition is preceded by large oscillations, while in some other runs, 
it occurs in one single step. 
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Fig. 5. A compilation of the variation of the normalized angular velocity of the 
liquid as a function of the imposed differential rotation rate Af. The velocity of 
the outer sphere is / ~ 4.5 Hz for all data shown. The shades give the proportion 
of data points that yield a given couple of values. Note that very different angular 
velocities are measured for a given Af. 

4-2 Evolution with the imposed differential rotation 



Figure [5] shows the dimensionless angular frequencies Af^ uid = U tp /27rsAf 
deduced from the electric potential differences measured at 40° latitude (ac- 
cording to equation [2]) , as a function of the imposed differential rotation Af, 
for a given ro tation rate f ^ 4.5 Hz of the outer sphere. The figure is similar 



to figure 4 of iNataf et all (120061 ). but data from several additional runs are 
included. The data points appear very scattered: for a given Af, the dimen- 
sionless angular velocity can vary by up to a factor of 3. We think that this is 
due to variations in the electric coupling between liquid sodium and the inner 
copper sphere. The coupling between the rotating inner sphere and the liquid 
sodium depends on the electric currents that are induced in the sodium and 
loop in the copper inner sphere. A similar scatter is observed when plotting 
the induced magnetic field as a function of Af (figure [6]). 

This scatter makes it difficult to draw robust conclusions about features such 
as super-rotation, but we will show in the next section that there is a linear 
relationship between the electric potential difference at 40° (AV400) and the 
azimuthal velocities in the sodium actually measured with ultrasonic Doppler 
velocimetry. We will therefore use AV400 as a proxy of the actual flow velocity, 
and we will see that several features of the flow become coherent when this 
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Fig. 6. A compilation of the variation of the radial component of the induced mag- 
netic field as a function of the imposed differential rotation rate Af. The velocity of 
the outer sphere is / ~ 4.5 Hz for all data shown. The colors indicate six different 
experimental runs. Note the large scatter of the data points, especially for negative 
A/. 
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Fig. 7. Angular velocity profile for / = 4.7 Hz and Af = 4.8 Hz. The time-interval 
between profiles is 30 ms, and there are 128 shots per profile. Probability density 
function or histogram of the velocity as a function of distance along shot. 

proxy is used in place of Af. 



4-3 Angular velocity profiles 



Using Doppler ultrasonic velocimetry, we are able to measure angular velocity 
of the liquid sodium flow. Figure [7] shows a profile of velocities measured along 
the shooting line for / = 4.7 Hz and Af = 4.8 Hz. Equation [3] shows that 
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Fig. 8. Velocity (in m/s) at a distance of 8.5 cm along the shot line (s ~ 13.6cm), 
as a function of A/ (left) and AV4o° (right). The rotation rate of the outer sphere 
is / ~ 4.5 Hz. The linear relation we get on the right shows that AT4o° is a good 
proxy of the angular velocity of the fluid at depth. 

the measured velocity component is proportional to the angular velocity of 
the fluid. Actually, the plot is built from a 40 s-long record and the shades 
map the probability density function (pdf) of velocity versus distance. In this 
regime, the flow is almost steady and the pdf is sharp. Velocity rises steadily 
from a value close to 0.3 m/ s at the outer boundary (top of the graph) to reach 
a plateau of about 0.85 m/s. The velocity further increases to a bump at 1.05 
m/s where the ultrasonic beam gets closest to the inner sphere (d ~ 18.8 cm, 
s ~ 7.9 cm). We interpret this bump as the signature of the magnetic wind, 
which is confined to the neighborhood of the inner sphere, where the magnetic 
field is largest. The plateau characterizes the region where the fluid is entrained 
at an almost uniform angular velocity. This profile is typical of those measured 
when the outer sphere is rotating, and differs from those obtained when the 
outer sphere is at rest. For all available profiles, we pick the velocity measured 
at a distance of 8.5 cm (on the plateau, and about at mid-depth of the fluid 
layer, in the equatorial plane, according to figure [T]) and plot it as a function 
of the corresponding Af and AV4o° (Figure [8]). While the former plot displays 
a large scatter, all points align well when velocity is plotted against AV400, 
illustrating the idea that the latter represents a good proxy of the effective 
velocity. From the regression line, we can translate AV4o° mTj ° fluid velocity 
or fluid flow frequency at depth in the fluid. We obtain: A/a U id/AV4o° — 2600 
Hz/V. Note that the fluid velocity deduced from AV400 using equation [2] agrees 
with the azimuthal velocity measured by ultrasonic Doppler velocimetry at the 
cylindrical radius corresponding to this latitude of 40°. 

In figure [9l we have plotted all the profiles for / ~ 4.5 Hz and positive A/, di- 
vided by the corresponding AV400. The horizontal scale has been chosen to pro- 
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vide the actual adimensional angular velocity for the reference case Af = 3.2 
Hz, according to equation [31 However, one should not consider this scale as 
universal, since the actual coupling is very variable (see figure E]). We have also 
drawn a synthetic profile derived from a numerical simulation similar to the 
one shown in figure [1] with Af = 0.4 Hz, but with the electric conductivity of 
the inner sphere reduced from its nominal value of 4.2 times that of the liquid 
sodium down to 0.04. In a first approximation, all experimental profiles col- 
lapse on a single curve, showing that the normalization with AV4o° is adequate. 
A closer look reveals that the rise of the velocity from the outer boundary is 
sharper when the forcing is weaker. In other words, the plateau of high angu- 
lar velocity extends further towards the outer sphere fo r weak forcing. We wi ll 
present in section l5\2l a simple model, in the spirit of iKleeorin et al.l (119971 ). 
that reproduces these characteristics (dotted and dashed lines). It seems that 
the relative size of the bump due to the magnetic wind above the plateau does 
not depend much upon the forcing. The comparison with the numerical profile 
suggests that the zone of strong magnetic wind is closer to the inner sphere 
in the experiments than in the simulation. We also note that the profile that 
corresponds to the highest differential rotation rate Af = 18.7 Hz plots below 
the others and has a distinct kink near the outer boundary, which makes it 
look more alike the profiles obtained with / = 0. 



4-4 Electric potentials 



We showed in iNataf et al.l (120061 ) that, when the outer sphere is at rest, the 
relative variation of electric potential differences (normalized to its value at 
40°) with latitude is independent of the differential rotation rate Af. The 
better data now available confirm this observation and further demonstrate 
that the latitudinal variation of potentials is close to, but significantly different 
from, the one expected for solid body rotation. Figure [TD] gives the variation of 
electric potential differences AV at latitudes of 10°, 20° and 30° as a function 
of AV400, for three different rotation rates of the outer sphere (/ ~ 4.5, 9 and 
13.5 Hz), and several different Af (both positive and negative). The case with 
no rotation (/ = 0) is included for reference. In all four cases, the data points 
align well, and fall not far from (but distinctly off) the straight lines calculated 
for solid body rotation (i.e., assuming U v proportional to s in equation [2]). In 
fact, we deduce from the measured trend for / = 4.5 Hz that the angular 
velocity at 30° is about 0.9 times that at 40°, in good agreement with the 
angular velocity profiles of figure M At lower latitudes, the apparent velocities 
deduced from the electric potentials are larger than the actual vel ocities, as the 
assumptions used in equation [2] break down (INataf et all 120061 ) . In contrast 
with this simple linear trend, there is a distinct interval of AV400 that displays a 
peculiar behavior, when the outer sphere is rotating. It occurs around AV400 = 
-1.7 mV for / = 4.5 Hz, -3.3 mV for / = 9 Hz, and -4.4 mV for / = 13.5 Hz. 
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Fig. 9. A compilation of the angular velocity profiles, normalized by AV4o°, f° r 
/ ~ 4.5 Hz and for imposed A/ ranging from 3.2 to 18.7 Hz. The horizontal scale 
has been chosen to give the actual adimensional angular velocity for the reference 
profile with A/ = 3.2 Hz. The continuous lines are predictions of different sort: the 
solid line is the profile obtained from the numerical model of figure (H (A/ = 0.4 
Hz) with the electrical conductivity of the inner sphere reduced by a factor 100 
(s ee section 15. 3p. The a lmost superposed dotted line is derived from the analysis 



of iKleeorin et al.l ()1997l ) , and the dashed line is the prediction of our extended 
torque-balance model for A/ = 15 Hz (see section [572]) . The two latter models only 
predict the shape of the boundary-layer profile, not the height of the plateau of 
near solid-body rotation. 

The kink means that, there, the variations in electric potentials are larger at 
low latitudes than at 40°. We will see below that the induced magnetic field 
also displays a very peculiar variation around the same values of AV400. 



4-5 Induced magnetic field 



The induced magnetic field is measured outside the sphere, in the labora- 
tory frame. The mean axisymmetric part of the induced field results from the 
interaction of the imposed axial dipolar magnetic field with the mean axisym- 
metric part of the meridional velocity field. The magnetometer signals are 
first averaged over 3 seconds, in order to remove t he oscillatory part (which 
is discussed in another article (ISchmitt et al.l . 120081 ) ) and other magnetic per- 
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Fig. 10. A composite of electric potential differences for four different rotation rates 
/ of the outer sphere. The electric potential difference measured at latitudes of 
10°, 20°, and 30° are plotted against that measured at 40°. The relation is almost 
linear, except for a kink observed around a negative value of AV4o°, which moves 
towards more negative values when / increases. The dashed straight lines indicate 
the relationship one would get if the fluid was in solid body rotation. Scattered 
points in the bottom panels are due to electromagnetic noise. 
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turbations. We then correct for temperature drift, by removing an empirical 
linear correction. In figure [HI both the radial B r and orthoradial Bg compo- 
nents are plotted as a function of AV4o° f° r several different imposed rotation 
frequencies / of the outer sphere. We have included the / = case for refer- 
ence. 

The first thing to note is the change of symmetry about AV400 = for small 
forcing : while the induced magnetic field is always positive when the outer 
sphere is at rest, we find that the induced field has the sign of AV400 (or A/) 
when the outer sphere is rotating. This can be directly related to the sense of 
the meridional circulation. When the outer sphere is at rest, the flow is always 
centrifugal in the equatorial plane, irrespective of the sign of Af. When the 
outer sphere is rotating, the sense of the mass flux through the outer viscous 
boundary layer is governed by the sign of the differential rotation between the 
fluid and the outer sphere (as for Ekman pumping). When Af/f is negative 
and small enough, the mass flux out of the boundary layer is thus centripetal 
at the equator. 



In that same regime, we also note that the induced magnetic field is weaker 
when the outer sphere is rotating : the meridional circulation violates the 
Proudman-Taylor constr aint and therefor e decr eases as the Ekman number 



decreases (see figure 6 of iHollerbach et al.l (120071 )). 



However, the most interesting feature of the curves shown in figure [TT] is the 
distinctive narrow peak observed at negative Af (here negative AV4o°). The 
peak occurs for a value of AV400 that moves toward more negative values as 
/ increases. In fact, the occurrence of the peak corresponds to the transition 
between two different regimes, which was invoked when discussing figure [TUJ 
We stress that it is because we plot the induced magnetic field as a function 
of AV4o° rather than Af that we obtain a well defined and reproducible peak 
(compare with figure[6]). The peak happens as the flow evolves from centripetal 
to centrifugal in the equatorial plane, as the effect of differential rotation takes 
over the effect of global rotation. 



4-6 Radial velocity profiles 



While the induced magnetic field measured outside the sphere is a good in- 
dicator of the global symmetries of the meridional flow, it cannot be used to 
investigate the internal organization of that flow. Here, we add the information 
brought by radial profiles of the radial velocity obtained by ultrasonic Doppler 
velocimetry. A compilation of the profiles obtained for three different rotation 
frequencies / of the outer sphere is given in figure [T2l for two different latitudes: 
—20° (figure [T2"r) and —40° (figure [T2"b). The profiles clearly confirm that the 
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Fig. 11. A composite of the induced magnetic field as a function of AV4o°, for four 
different rotation frequencies / of the outer sphere. Radial component on the left, 
and orthoradial component on the right. Note the symmetry change around the 
origin (A/ ~ 0) between the case without rotation (/ = 0) and the other cases. A 
striking peak is observed for a particular value of A/ in the counterrotating regime 
(AV400 < 0). The colors correspond to different runs. 



sign of the radial velocity changes with the sign of A/, as inferred from the 
induced magnetic field measurements. At —40°, the profiles are characterized 
by a broad peak with radial velocity rising from zero at the inner sphere to a 
maximum, which can reach about 3% of the tangential velocity of the inner 
sphere. There, the flow is centripetal when the Rossby number is positive. The 
radial velocity falls back to zero at nearly the same radius for all /, but the 
profiles for the higher / (lower Ekman number) are characterized by a strong 
second peak. At —20°, the amplitudes are similar but seem to increase with 
/. The change of sign as Af changes sign is well visible again. 



22 




adimensional radial velocity at -20° adimensional radial velocity at -40° 



Fig. 12. Radial velocity, obtained by ultrasonic Doppler velocimetry, as a function 
of adimensional radius (vertical axis) at two different latitudes: —20° (a) and —40° 
(b). The dotted line indicates the radius of the tangent cylinder on the inner sphere, 
whose adimensional radius is 0.35 (bottom of the plots). All velocities are normalized 
with the differential tangential velocity at the equator of the inner sphere. The 
profiles are labelled by their (/, A/) pair. Note that the sign of the radial velocity 
changes as the sign of Af, as expected for Ekman-pumping dominated flow. 



5 Discussion and conclusions 



5. 1 Moderate-Rossby number regime 



The magnetostrophic regime at moderate Rossby numbers is characterized by 
a low level of turbulence. Most of the energy is in the axisymmetric mean 
flow described in this article. The meridional circulation has the symmetry 
of Ekman pumping: counterclockwise in a meridional plane (i.e., centrifugal 
at the equator) for Ro > and clockwise for Ro < 0. Radial velocities are 
at least one order of magnitude smaller than azimuthal veloci ties. Waves are 



obser ved in this regime, and are described in another article ( iSchmitt et al. 
2008h . 



Profiles of the angular velocity of the flow have been obtained, using ultrasonic 
Doppler velocimetry. They reveal a region where magnetic wind is present, 
near the inner sphere, where the magnetic field is largest. Farther away, the 
flow is geostrophic, under the effect of the Coriolis force. A large part of the 
geostrophic region is entrained by the magnetic wind at a nearly constant 
angular velocity. 
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The measurements show that the shape of the angular velocity profile is nearly 
independent of Ro as long as it remains moderate. This suggests that, in con 



trast to the case with a fixed outer sphere (INataf et all 120061 ; iHollerbach et al. 
20071 ). the linear solution (Ro << 1) provides a good description of the flow. 



We were able to compute solutions in the linear regime with the geometry 
and parameters of the DTS experiment for / = 5 Hz (E = 4.7 10~ 7 ), as 
shown in figure HJ In figure [1^, one clearly sees a non-geostrophic inner region 
where the Ferraro law of isorotation applies, a central part where the flow is 
geostrophic and in almost solid-body rotation, and an outer domain, which is 
also geostrophic but in which the angular velocity decreases rapidly towards 
the outer sphere. In order to compare with the experimental results, synthetic 
velocity profiles have been computed. The comparison of the angular velocity 
profiles is shown in figure [9j All profiles show the plateau of nearly constant 
angular velocity, the central peak due to the magnetic wind, and the gentle 
drop to lower velocity towards the outer sphere. We note however that this ve- 
locity drop is much steeper in the numerical simul ation (which ha s Ro = 0.08). 
Several of these features had been predicted by iKleeorin et al.l (119971 ) in an 
asymptotic study of spherical Couette flow in a dipolar magnetic field. We 
recall their main results below and propose an extension of their approach, 
which better explains our experimental results. 



5.2 Modified Taylor constraint and extended quasi-geostrophic model for the 
moderate-Rossby number flow 



The quasi-geostr ophic regime wi t h an imposed magnetic field has been stud- 
ie d extensively by IKleeorin et al.l (119971 ). Following the pionneer computations 
of IHollerbach! (119941 ) of rotating spherical Couette in a dipolar magnetic field, 
they worked out the asymptotic solutions for small Elsasser, Ekman, Reynolds 
and Rossby numbers. They predicted the existence of a geostrophic region 
where the magnetic torque on coaxial cylindrical tubes balances weak fric- 
tion in the upper and lower Ek man layers, thu s providing an extension of the 
zero-Ekman number analysis of iTaylorl (119631 ). which results in the so-called 
Taylor-state. We consider the torques at work on a cylindrical tube at cylin- 
drical radius from s to s + ds. Neglecting the volumetric viscous friction (small 
Ekman number), the torque balance can be written as: 



d^ T (s) + d^ B (s) + dT M (s) = 0, 



(6) 



where djT and d'jB are the (resistive) torques on the top and bottom spherical 
caps at the outer sphere, and oTm is the volumetric (driving) torque exerted 
by the Lorentz force (magnetic torque). In our case, the surface torques are 
due to Ekman friction on the outer sphere and can be written: 
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dj T (s) = dj B (s) = 2irs 2 T V -, (7) 

COS u 

where 9 is the colatitude and ry is the azimuthal component of the friction 
stress on the outer shell. 

At high latitudes, friction in the Ekman layers is weak and, as in the Taylor 
state, the balancing magnetic torque vanishes in the steady-state and the fluid 
is in solid body rotation. The magnetic torque can be expressed as: 



r ZT -> -> 2,71 d 

dTM(s) = 2iis 2 ds / (j x B) ip dz = — ds— 
Jzb j^o ds 



B s b v dz 



where B s is the component of the imposed magnetic field along the cylindri- 
cal radius s, while b v is the (unknown) azimuthal component of the induced 
magnetic field. This component results from the shearing of the dipolar field 
lines by deviations of the fluid flow from solid body rotation. A weak mag- 
netic torque thus implies a small deviation from solid body rotation, especially 
where the imposed field B s is large, i.e. near the inner sphere (b ut not too close, 
becau se the flow is no longer geostrophic there) . As noted by iKleeorin et al 
( 119971 ). this explains why the angular velocity of the fluid Au(s) is almost 
constant in the central part of the spherical shell. The progr essive decrease 
of the angular velocity as s increases further was predicted by IKleeorin et al 
(119971 ). In this layer, which they call "magnetic-Proudman" , shear in the up- 



per and lower Ekman layers increases, so as to produce a large enough b v to 
compensate for the decrease of B s in the expression of the magnetic torque. 
In order t o be more quantitatiy e, one needs to express b v (s,z) as a function 



of Au(s). IKleeorin et al.l (119971 ) show that, when the imposed magnetic field 



is a dipole, the magnetic torque can be written: 



d T M ( s ) = l^RVAft- 



2 1 



3/2 Au ] _ A ^( g ) 



(9) 



where Ao;' is the differential rotation of the fluid in the region of nearly solid- 
body rotation. The friction stress for the laminar Ekman layer is easily calcu- 
lated and yields: 



d^ris) = — 2irpa A AVLv vVL 

a a 6 



1-1/4 



Auj(s) 



(10) 



Injecting express ions 191 and [TUl i n equ ation [6] in the boundary layer approxi- 



mation (s ~ 
angular velocity versus s: 



Kleeorin et al.l (119971 ) obtain the shape of the profile of the 
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The synthetic angular velocity profile derived from this expression is the dot- 
ted line of figure [9j The value of Au^ has bee n chosen to coincide with the 
full numerical profile drawn (solid line). Indeed. iKleeorin et al.l (119971 ) studied 
the case when the flow remains geostrophic all the way to the inner sphere, 
where they matched the solid-body solution to the Stewartson layer modified 
by the magnetic field. They could thus predict the value of AoA In our case, 
the inner region is not geostrophic, and we only use their analysis to probe 
the shape of the angular velocity profiles in the ge ostrophic region . It is clear 
that, in that region, the asymptotic prediction of IKleeorin et al.l (119971 ) is a 
very good approximation of the actual numerical profile. However, we note 
that both profiles fail to explain the much more gentle rise of angular velocity 
observed in the experiments. The explanation is that the Ekman boundary 
layer becomes turbulen t at the high A f reac hed in our experiments. We can 
retain the approach of IKleeorin et al.l (119971 ) and simply allow for turbulent 
friction stress ry in equation [7J Because of turbulence, Ekman friction is en- 
hanced and a stronger magnetic torque is needed to balance it. This results 
in a widening of the "Proudman-magnetic" layer, as shown by the dashed 
line in figure O computed for Af = 15 Hz. The widening increases as Af 
increases, in excellent agreement with our measurements. The full ingredients 
and implications of this extended model will be given elsewhere. 



5.3 Variable electric coupling and Ro e ff 



In iNataf et al.l (120061 ). there were hints that different flows could be obtained 
for the same imposed rotation rates / and Af. Our new data confirm this 
trend. We speculate that this behavior is due to variations in the electric cou- 
pling between liquid sodium and the copper inner sphere, under the effect of 
electrowetting feedback. It is thus difficult to relate the observed signals to 
the imposed Af, as we would like to do in order to study the possibility for 



the m agnetic wind to induce super-rotation of the fluid as in iDormy et al. 



(119981 ). However, we find that most results can be cast in a coherent frame if 
we use AV4o°, the electric potential difference at latitude 40°, as a proxy of 
the actual angular velocity, in place of Af, resulting in an effective Rossby 
number Ko e ff. In order to mimic a weaker electric coupling, numerical sim- 
ulations have been run where the electric conductivity of the inner sphere is 
reduced as compared to that of copper (nominal conductivity ratio of 4.2). 
Synthetic angular velocity profiles are shown in figure [T3k for a conductivity 
ratio of 0.04 and 0.01. The profiles retain the same shape but the overall am- 
plitude is severely reduced. This explains why we could get coherent results by 
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Fig. 13. Synthetic Doppler velocity profiles computed from numerical simulations, 
(a) angular velocity synthetics for a shot from the azimuthal assembly at 10° (see 
projected raypath in figure [Ih.). The synthetics take into account the finite width 
of the ultrasonic beam and the projection of the meridional velocity. The dash-dot 
line is for the numerical simulation of figure [U in which the electric conductivity 
of the inner sphere is 4.2 times that of the liquid sodium, as expected for copper. 
The other computations have a reduced inner sphere conductivity ratio of 0.01 
(dot-dot-dash), and 0.04 (solid line). The corresponding angular velocity profile for 
the latter is also drawn in figured! (b) radial velocity synthetics for radial shots at 
latitudes ±40°, ±20° and ±10° for the latter model. Velocities are normalised with 
the differential tangential velocity at the equator of the inner sphere, and distances 
are normalised by the radius of the outer sphere. 



renormalizing all data with the AV4o° proxy. Comparing the spread in figure 
[5] with the numerical simulations suggests that the effective conductivity at 
the sodium/copper interface could be reduced by as much as 100 by wetting 
problems. We observe that the electric coupling deteriorates with time over 
the months, but that it can be somewhat revived by letting the copper in 
contact with liquid sodium at 150°C during 3 days. 

The numerical simulations show that the meridional circulation is not very 
sensitive to the conductivity ratio. Radial profiles of the radial velocity are 
shown in figure fl~3b for the numerical model with the conductivity ratio of 
0.04. As in figure [T)d, the circulation is centrifugal at the equator (Ro > 0) 
and the radial velocity is positive at the latitude of 10°. It is negative at —40°, 
in good agreement with our measurements (figure fT2b). At 20°, the profile is 
more complex, but it is also negative near the inner sphere, as experimentally 
observed (figure [T2h). In the numerical model, the radial velocities reach about 
1% of the inner sphere tangential velocity, while it is about twice as large in the 
experiments. The experiments show clear evidence of a secondary meridional 
cell farther away from the inner sphere. 
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5.4 The peculiar regime for Ro e jf ~ — 1 



In the counter-rotating case, the centripetal circulation at the equator, which 
characterizes the moderate Rossby number regime, transforms into a centrifu- 
gal circulation at the equator as the Rossby number increases. We recall that 
at small Rossby number, the meridional circulation is governed by Hartmann- 
pumping at the inner sphere, and is therefore directed towards the inner sphere 
when it is counter-rotating, while centrifugation from the inner sphere takes 
over at large Rossby numbers. The transition between the two regimes occurs 
over a narrow interval of AV4o° > and is marked by a peak of the induced mag- 
netic field. This happens where the effective angular velocity of the fluid is 
such that it almost cancels that of the outer sphere, as seen in the laboratory 
frame. 

Our interpretation is that, in a narrow window of AV400, the antagonistic 
effects of the rotation of the outer and inner spheres drive the liquid in a 
state that shows almost no rotation with respect to the laboratory frame. The 
constraint of rotation then vanishes, and a vigorous meridional circulation is 
allowed, resulting in an enhanced induced magnetic field. 

Some support for this interpretation is given by a direct measurement of the 
angular velocity using the ultrasonic Doppler technique. As shown in figure dU 
a peculiar profile is measured close to the B peak. An almost uniform angular 
velocity is measured (down to d ~ 7.4 cm, i.e. at mid-depth of the fluid shell). 
Using equation [31 we find a fluid rotation frequency of -5.7 Hz on this plateau. 
Since the ultrasonic probe is installed in the rotating sphere, the angular 
velocity with respect to the laboratory frame is obtained by adding the angular 
velocity of the outer sphere. The resulting angular velocity is then -5.7 + 4.2 
= -1.5 Hz, which is small but non zero. Note that the electric potentials (figure 
[TOT) also suggest that, on and beyond the peak, the fluid is nearly in solid-body 
rotation. However, the above scenario is probably too simple, and the variation 
of the induced magnetic field calls for a more sophisticated interpretation. For 
example, the B r and Bg records of figure [3] show that the "centrifugal" regime 
is announced by a strong rise of Be, while B r starts decreasing instead. Besides, 
pure solid-body rotation would imply no meridional circulation altogether. 



5. 5 Implications for experimental dynamos 



The transition between the centripetal and the centrifugal regimes deserves a 
sp ecific discussion in th e context of the preparation of experimental dynamos. 



In ICardin et al.l (120021 ). we evaluated the possibility to run a dynamo where 



the forcing is a rotating spherical Couette flow. Such a dynamo would run in 
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Fig. 14. A profile (pdf) of the angular velocity of the fluid for / = 4.2Hz and 
A/ = — 18.2Hz, for which AV40 = — l-7mV. This value corresponds to the peak 
observed for the induced magnetic field. 



the magnetostrophic regime, expected for planetary cores, and it would be of 
the aw-type, displaying large azimuthal magnetic fields. However, the crit- 
i cal magnetic Reynolds num ber Rm c for these dynamos appears to be high 
(jSchaeffer and Cardinl . 120061 ). Indeed, there is an opti mum ratio of merid- 



ional to azimuthal v elocities that yield a minimum Rm c (IRavelet et all 12005 



Bayliss et all 120071 ). Here, we find that, in a rapidly rotating spherical Cou- 



ette flow, the meridional velocities are at least one order of magnitude smaller 
than the azimuthal v elocities. In contrast , the mean flow in ot her experiments 



(IPefflev et all 1200(1 IRavelet et all 12003: iBavliss et all l2007h is tuned to be 
close to an axi-symmetric flow of the Dudley and Jamesl ( 1989 ) family, with 
a strong meridional c i rculation, yielding an q? 2 -dynam o, and moderate Rm c 



(IKumar and Roberta . 1 19751 ; Dudley and Jamesl . Il989l ). It has been recently 
suggested that the dynamo efficiency of such flo ws can be significantly reduced 



by turbulent fluctuations (IBavliss et all 120071 ). Now, in the narrow window 



that corresponds to the transition from the centripetal regime to the centrifu- 
gal regime, the magnitude of the meridional circulation becomes comparable 
to the azimuthal one. Thus, the value of the magnetic Reynolds number for 
the onset of dynamo action may be minimum around this transition. A rotat- 
ing spherical Couette flow dynamo experiment should thus make it possible 
to observe the transition from an a 2 -type (with Rm c perhaps large compared 
to its laminar value) to an aw-type dynamo (with a large Rm c for the lam- 
inar flow but weak turbulent effects). We are currently conducting extensive 
numerical simulations of the dynamo instability for spherical Couette flow to 
document this effect. 
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5. 6 Implications for core dynamics 



The modified Taylor state observed in the DTS experiment bears some rel- 
evance to the expected dynamical regime in the Earth's core. Indeed, the 
Ekman number based on the expected fluid viscosity is very small, but obser- 
vations of the nutations reveal that the friction at the core -mantle boundary is 



some 10 larger t han expected for linear Ekman friction (jHerring et all 12002 



Koot et al. . 2008h. Although t he origin of this enhanced fricti on is still debated 



(IDeleplace and Cardinl . 120061 ; iBuffett and Christensenl . 120071 ) and is undoubt- 



edly different from the turbulent friction present in DTS, it would be worth in- 
vestigating its effect on the dynamics of the core, fr om a modified Taylor state 
point of view. Quasi-geostrophic numerical models (jSchaeffer and Cardinl . 120061 ). 
in which Ekman friction is parameterized, are a particularly attractive tool 
for quantifying this effect. 
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